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Abstract 

High precision spectroscopy can provide a sensitive tool to test Coulomb's law on atomic 
length scales. This can then be used to constrain particles such as extra "hidden" photons 
or minicharged particles that are predicted in many extensions of the standard model, 
and which cause small deviations from Coulomb's law. In this paper we use a variety 
of transitions in atomic hydrogen, hydrogenic ions, and exotic atoms to probe Coulomb's 
law. This extends the region of pure Coulomb's law tests to larger masses. For hidden 
photons and minicharged particles this region is already tested by other astrophysical and 
laboratory probes. However, future tests of true muonium and muonic atoms are likely to 
probe new parameter space and therefore have good discovery potential for new physics. 
Finally, we investigate whether the discrepancy between the theoretical calculation of the 
^ s i/2 1 ~ ^l/2 2 transition in muonic hydrogen and its recent experimental measurement at 
PSI can be explained by the existence of a hidden photon. This explanation is ruled out by 
measurements of the Lamb shift in ordinary hydrogen. 

1 Introduction 

Precision spectroscopy has a long standing record of providing insights into fundamental physics. 
In particular, the discrete nature of spectral lines has led to the discovery of quantum mechanics 
and the discovery of the Lamb shift was one of the first confirmations of quantum electrody- 
namics. With its ever increasing precision, spectroscopy can continue to provide a powerful 
probe of new physics. 

Concretely, in this paper we want to use spectroscopy to test Coulomb's law with high 
precision on atomic length scales [THS]. This in turn allows us to obtain constraints on new 
particles such as hidden photons [HH1E] or minicharged particles [Q[7] which arise naturally in 
a variety of extensions of the standard model [H 151-113] (see also [13] for a review £]■ 

At this point let us note that tests of Coulomb's law are an especially clean and model 
independent probe of such new particles as the sensitivity does not depend on the stability of 
the particles or the presence/absence of certain decay channels. To illustrate this point, let 
us take fixed target experiments (see, e.g. [16]) as an example. These experiments typically 
rely on a displaced vertex and therefore on a relatively long decay length of a hidden photon 
produced in electron-nucleus collision. Now, let us assume that there is also matter charged 
under the hidden U(l) (i.e., it couples directly to the field in Eq. (|2.ip ). If this additional 
"hidden matter" is lighter than my /2, the hidden photon decays quickly into a pair of these 

Spectroscopy can even be useful for constraining Unparticles (see e.g. [15]). 
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particles, dramatically shortening the decay length and possibly invalidating the bound. Tests 
of Coulomb's law therefore provide interesting complementary information. 

In the following our prime example will be hidden photons which (as we will briefly recall 
below) cause a deviation of Coulomb's law of the form 

V(r) = - — (l+e~ m -r'V) (1.1) 
r 

where my is the mass of the hidden photon and x is a so-called kinetic mixing [8]. Independent 
of this particle interpretation, our bounds can also be taken as a parametrisation of a deviation 
from Coulomb's law by a Yukawa type potential with a characteristic length scale ~ 1/my. 
The length scale of ordinary atoms are typically of the order of the Bohr radius and correspond 
to masses of the order of keV. Exotic atoms, such as muonic hydrogen where the electron is 
replaced by a muon are even smaller, being most sensitive to masses of order MeV or above. 
To complement existing spectroscopic tests [IH3] we will pay particular attention to the latter 
regime. 

The paper is set up as follows. In the next section we will briefly recall how new particles 
can modify Coulomb's law focusing on our main example, the hidden photon. Minicharged 
particles are discussed in App.[A] In section [3] we give a detailed explanation of the method we 
use to constrain hidden photons, including a discussion of how to obtain the proper behaviour 
for small (and large) hidden photon masses (section 13. 1 . 2|) . In section [3] we apply this method 
to construct bounds from a variety of atomic spectra. Sections 14.11 and 14.21 deal with ordinary 
hydrogenic atoms, whereas sections I4.3U4.4| and 14.51 deal with exotic atoms. Particularly inter- 
esting speculative bounds are calculated for true muonium and muonic hydrogen. In sections 
14.4.21 and 14.51 we show that future measurements of true muonium and muonic atoms promise 
significant discovery potential for new physics. Finally, in section [5] we discuss whether the 
hidden photon can be used to explain an anomaly recently observed in muonic hydrogen. We 
conclude in section [6l 



2 Hidden photons and deviations from Coulomb's law 

Hidden photons can interact with the photon via a so-called kinetic mixing [HE], 

L = -\f^ - + |iV^ + IrrfyBuB" (2.1) 

where F^ u and are the photon (A u ) and hidden photon (B u ) field strengths. Both 
and B^y are individually gauge invariant, so the third mixing term is allowed. Since the mixing 
term has a mass dimension of 4, the kinetic mixing parameter x nas a mass dimension of zero 
and therefore is a renormalisable coupling. This means that it is not suppressed by any higher 
mass scales and should be observable (even if small) at lower energies. We also note the fourth 
term which indicates that the hidden photon has a non zero mass. From the point of view of 
low energy effective theory, x an d my are simply free parameters which are not constrained by 
any particular physical mechanism. However, extensions of the standard model based on, for 
example, string theory predict values in the range [9lfT2|fT3]. 

10~ 12 < X < HT 3 . (2.2) 

A summary of the current bounds on the parameter x can be found in Fig. [TJ 

The kinetic mixing term in Eq. (|2.ip causes a tree level insertion to the photon propagator as 
shown in Fig. [2J An important physical implication of this is the addition of a new Yukawa-type 
term to the Coulomb potential, 

V(r) = - — (1 + e""V V) = Fcouiomb(r) + 6V(r), (2.3) 
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Figure 1: Summary of the current bounds on hidden photons (compilation from |14pi7| updated 
with OS]). We note that in addition to the bound labelled "Coulomb", the bounds labelled 
"Earth", "Jupiter" and "Rydberg" also originate from tests of Coulomb's law. The best bounds 
that we derive from atomic spectra are represented by the black lines. The dashed black line is 
from the Lamb shift in atomic hydrogen. The dotted black line is a combined bound from the 
lsi/2 — 2si/2 an d 2si/2 — 8si/2 transitions in atomic hydrogen. The solid black line is the bound 
obtained from the Lamb shift in hydrogen-like helium ions. The light orange areas correspond 
to regions suggested by astrophysical and cosmological puzzles |18H22| . The brown region is 
derived from measurements of the Rydberg constant, and represents bounds already obtained 
from atomic spectra [TH3]. 




Figure 2: Feynman diagrams for the potential between two charged particles, Eq. (|2.3p . In 
particular the second diagram gives the hidden photon contribution to the interaction between 
two charged particles, leading to a modification of Coulomb's law. 
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where Z is the charge of the massive central particle which causes the potential and a is the 
fine structure constant. 

Note the following limits, 

• my — > 0: 

Here the exponential term tends to unity and we recover the original Coulomb potential 
up to a factor (1 + X 2 )- This can be absorbed in the definition of a, making it unobserv- 
able. Physically this is sensible, since in the zero mass limit our hidden photon becomes 
indistinguishable from the standard model photon, at least as far as the electromagnetic 
force is concerned. Therefore no physical effects, such as a deviation from Coulomb's law, 
should be observable. 

• my — > oo: 

Here the second exponential term dies off and leaves us with the original Coulomb po- 
tential. Again this makes sense, since in the large mass limit the hidden photon becomes 
impossible to excite as a virtual particle and therefore should not contribute to the elec- 
tromagnetic force. 

• Intermediate my: 

We expect non zero deviations from Coulomb's law in this region. In particular, we will 
find that the most significant deviation occurs for my ~ ^ where Iq is the length scale of 
the relevant physical process. 

Therefore any constraints on hidden photons based on deviations from Coulomb's law will 
be strongest around my ~ ^ and drop off on either side in the low and high mass limits. Here, 
Iq will be given by the length scale of the atom in question. There exist a range of different 
atomic systems, including ordinary hydrogen, muonic atoms, and more exotic atoms. Therefore 
we can test a wide range of masses, from keV to larger than MeV. 

3 Obtaining constraints using atomic spectra 

To obtain our constraints we adapt the method presented in Ref. [6], where the measurement 
of the Lamb shift is used to derive a bound on minicharged particles. 

At first order in perturbation theory the energy shift of a state \tp n ) is given by 

5E^ = (rp n | R> | = Q, n | SV | Vn> (3.1) 

where the \ip n ) are taken to be the 0th order wave functions. For this to be a good approximation 
the energy shift should be small. This is consistent with what we expect, since so far no large 
deviations from the standard QED predictions have been observed. If the standard prediction 
and the experimentally measured values agree, then we can impose that SEn^ must be smaller 
than the uncertainty in the transition. This will result in a bound on 5V. 

Let us briefly comment on some points relating to these uncertainties that will be relevant 
for our discussion. 

• "Same n" and "different n" transitions 

We can write the theoretical energy of an arbitrary state as; 

E n ,l,j = En,] R + ^n,l,j (3-2) 

where the first term is the sum of energies from the Dirac equation plus recoil corrections 
(effectively the 0th order energy). The second term is the Lamb shift (defined as any 
contribution which separates states of the same n,j). 
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The first term is proportional to the Rydberg constant i?oo and therefore will have an 
uncertainty of approximately 10 -10 eV [23]. This means that any transitions between 
states of different n in atomic hydrogen will have a theoretical error contribution of around 
10" 10 eV from R^. 

The situation is worse in exotic atoms, as the effective Rydberg constant is modified by a 
factor proportional to the reduced masses ^exotic of atomic hydrogen and the exotic 
atom respectively, 

D /^exotic j~f fo o\ 

K e ff — /too. (0.6) 

So that in, for example, muonic hydrogen, there would be a fractional uncertainty of 
around 10 due to the mass of the muon, which would cause an overall uncertainty of at 
least 10" 4 eV. 

Note that transitions between same n states do not have a 0th order energy and are 
limited only by uncertainties in the Lamb shifts of the states. 

• Definition of "uncertainty" 



To be very conservative when forming bounds, we estimate the total uncertainty of a 
transition by adding together the absolute values of the theoretical and experimental 
errors, i.e. for a given measurement M we use 

AM = \AM(th)\ + \AM(exp)\ (3.4) 



The theoretical contribution to AM will come mainly from uncertainties in the finite 
nuclear size, with an additional contribution from the Rydberg constant in the case of 
different n transitions. For most of the atoms we consider, data for the finite nuclear size is 
coherent. For example the hydrogenlike helium ion has three different electron scattering 
determinations, and also a muonic helium ion determination, of the alpha particle charge 
radius. These values all agree within 1 a |24j . meaning that to a good approximation we 
can calculate the theoretical value for a transition by assuming one particular value of the 
nuclear radius. We can estimate the uncertainty from finite nuclear size effects simply by 
considering the uncertainty in this one value of the nuclear radius. This is indeed what is 
done in theoretical calculations to which we refer, and the corresponding error is included 
in \AM th \. 

However when we consider atoms with a proton nucleus things become more compli- 
cated. The recent muonic hydrogen determination of r p = 0.84184(67) fm [25] gives us 
the most precise measurement from atomic spectra. This disagrees with our best previous 
atomic spectra determination of r p = 0.8768(69) fm from [22] by around 5 a. The muonic 
hydrogen extraction also deviates by around 2.5 a from the best electron scattering de- 
termination of r p = 0.897(18) fm [26]. To be conservative we therefore modify our error 
analysis to take into account the large variation in r p . We do this by adding an additional 
term \AM(r p )\, accounting for the discrepancies in the r p measurements, 

AM* = \ AM{th)\ + |AM(ea;p)| + \AM(r p )\. (3.5) 

To leading order the proton radius contribution to a given state is 



4^2 
V 



ENs(r p ) = — ^3 — S m , (3.6) 

where m is the mass of the orbiting particle. In atomic hydrogen m = m e , and in 
muonic hydrogen m = m«. We make a rough estimate for \AM(r p )\ by considering 
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two widely separated values of the proton radius. Denoting r Pt n = 0.84184(67) fm as 
the muonic hydrogen determination, and r p ^ e = 0.897(18) fm as the electron scattering 
determination, we set 

\AM{r p )\ = \E NS (r p , e ) - E NS (r p>li )\. (3.7) 

To form bounds we use AM* for atoms with a proton nucleus, which in our paper are 
atomic hydrogen and muonic hydrogen. We use AM for all other atoms. 

Note that errors quoted in the text will always be at the 1 a (~ 68% confidence level). 
However, unless otherwise stated we will consider 2a errors (~ 95% confidence level) for 
the bounds produced in the figures. This essentially means that all terms on the right 
hand side of (|3.5|) . (|3.4j) are multiplied by a factor of 2. 

• Radius of the proton and other nuclei 

As already discussed the (charge) radii of the nuclei are a major source of uncertainty. In 
addition, we have to take care that our determination of the radii is from an independent 
source. For example we can not take a measurement of the Lamb shift in ordinary 
hydrogen to measure both the radius of the proton and put a bound on deviations from 
Coulomb's law. Two independent measurements are needed. Moreover to avoid even 
partial degeneracies (which tend to weaken the bound in particular at short length scales) , 
it is best if the determination of the radius is obtained at relatively high momentum 
transfer, corresponding to a short length scale. For this reason we will mainly use the 
values obtained from electron scattering data. 

3.1 Searching for deviations from Coulomb's law 
3.1.1 Naive bounds 

Let us briefly demonstrate how this method works by using the Lamb shift, i.e. the 2s-yji — 2pi/2 
transition in atomic hydrogen to constrain the hidden photon. We use 5V given by Eq. (|2.3I) , 
For this case Eq. (|3.ip evaluates, 

SE = J drr 2 5V(r) [R 2 20 (r) - R 2 21 (r)] = - X 2___Z_ (3 . 8) 

where we have used the normalised radial hydrogen wave functions R n £ given by 

Rw{r) = 71 ^ v 1 " f ) exp ^ 2 ) ' R ^ = ^/E^h p ex p(-^/ 2 ) < 3 - 9 ) 

and where p = r/a and a -1 = am e . 

We use an experimental uncertainty of 3 kHz [27] and theoretical uncertainty of 6 kHz |28j , 
as well as a contribution of 17 kHz from \AM(r p )\ ([377]) ■ We use ([33]) to get AM* = 10~ 10 eV 
at the 1 a level. Note that the (blue) curve for the 2s]/ 2 — 2pi/ 2 transition (Fig. [3]) is at the 2 
a level. 

This has the correct shape; the bound dies off the in the limits my — > and my — > oo, and 
is strongest at my ~ -, where the Bohr radius a is the typical length scale involved. 

We can do the same for the lsi/2 — 2si/2 transition in atomic hydrogen. The experimental 
value has a relative uncertainty 2.8 x 10 -14 and represents the most precise measurement of 
atomic hydrogen [29j. However, the bounds are limited by a much larger theoretical uncertainty. 
As already discussed above there are uncertainties of around 10~ 10 eV from the Rydberg con- 
stant as well as a similar contribution from the Lamb shifts of the states [30J. We also need to 
add the |AM(r p )| contribution of ~ 7 x 10~ 10 eV to get AM* = 1 x 10~ 9 eV at the 1 a level. 
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Figure 3: The blue curve denotes the bound on hidden photons obtained from the 2s]/2 — 2>Pi/2 
transition in atomic hydrogen. We use a conservative error of AM* = 2 x 10~ 10 eV at the 2 
a level where AM* is defined in (|3.5p . Note that for same n transitions x dies off correctly 
for both small and large my. The red curve shows the naive bound from the lsi/2 — % s i/2 
transition in atomic hydrogen. It has an incorrect behaviour at small my. The green curve 
gives the correct bound obtained by combining the lsi/2 — 2si/ 2 with the 2si/ 2 — 8S1/2 transition, 
according to the procedure described in Sect. 13.1.21 The green and blue bounds will turn out 
to be the best ones that we can derive from atomic spectra, and are combined together in Fig. 
[TJ For comparison we depict earlier bounds on hidden photons as the colour filled regions. 
Those corresponding to pure tests of Coulomb's law (also obtained from atomic spectra) are 
highlighted in brown [T]-[3]. The remaining white region corresponds to unexplored parameter 
space. 



For the first order energy shift we have 

roc 

5E = drr 2 5V{r)[Rl Q {r) - R\ Q {r)] (3.10) 

J o 

12 + amy (60 + amy (87 + 14amy (4 + amy))) 

4a(l + amy) 4 (2 + amy) 2 \ ' 

where 

2 

R ™( r ) = ^3/2 ex P(~^)- ( 3 - U ) 

The result is shown in Fig. [3] as the red curve. This bound does not have the correct drop 
off for small masses. We can understand why by looking at SV(r) in Eq. (|2.3p . which dies 
off at large masses but grows at small masses. Therefore we get a bound which saturates at 
small masses, which is not physically correct. This is simply an artefact of the splitting of 
the potential that we have chosen to set up our perturbation theory. At small masses our 
perturbation reduces to a term that has the form of a Coulomb potential, but with an extra 
factor (1 + x 2 )- This effectively increases the strength of the electromagnetic coupling and 
therefore the energy difference between the two states. In other words we have forgotten to 
properly (re-)normalise the coupling a, i.e. we have not absorbed the factor (1 + x 2 ) into our 
definition of a. In the following we will show explicitly how this can be remedied. 

One might wonder why this problem does not affect the Lamb shift bound. The reason 
is simply that for a perfect Coulomb's law the energies of 2si/ 2 and the 2pi/ 2 are degenerate. 
Therefore, adding a term of the Coulomb's law form does not produce an energy shift between 
the two states at lowest order. This is true for all same n transitions. 



= X 2 a 
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3.1.2 Bounds including a proper renormalisation of a 

As we have seen above, the problem is that for both my — > oo and my — > the potential 
Eq. (|2.3|) has exactly the 1/r behaviour of a Coulomb potential. However the coupling constant 
differs by a factor (1 + x 2 )- a becomes a function of x an d my, i.e the hidden photon alters 
the coupling constant itself. In that sense a becomes an unknown quantity that needs to be 
fixed by experiment. Since we now have two unknowns, x an d ct, we need two measurements to 
solve for them (of course the same strategy also works if we allow other / additional parameters 
to vary). We will briefly sketch how we do this. 

Let us assume we have two observables M\ and M 2 . Theoretically these are functions of a, 
X 2 and my. To keep the notation transparent we will suppress the dependence on my in the 
following. Therefore, we have Mi(a,x 2 ) and M 2 (a,x 2 )- 

Now we have two measurements and results are often quoted in the form 

M l I cxp - Mi|th = SMi ± AMi, (3.12) 
M 2 | eX p - M 2 |th = $M 2 ± AM 2 , (3.13) 

without considering a hidden photon. Therefore, in our setup this means 

MxUp - Mi(a , 0) = SMi ± AM 1; (3.14) 
M 2 | cxp - M 2 (a , 0) = 6M 2 ± AM 2 , (3.15) 

with some value ao such that both 5M\ and 5M 2 are small. If ao can be chosen that 8M\ and 
5M 2 are compatible with within the errors, then the measurements are consistent with the 
standard model and no hidden photon is required. 

Now, what happens if we include a hidden photon? As we are interested in small corrections, 
we can expand about (a, x) = (c*o, 0). Hence we obtain 

^^-\ a=aoa 2 =Q 6a + ^^-\ a=ao , x 2=oX 2 = SMi ± AMi, (3.16) 
dM 2 dM 2 2 

-7^-\ a=aOiX 2 =0 6a + -^2-| Q=Q0iX 2 =0 x = SM 2 ± AM 2 . (3.17) 
In matrix notation this linear system of equations reads, 



dM_i\ 8Mi 

\a=a ,x z = 
-Q^r\a=a ,x 2 =0 



^' - - 9 -hL=^=o U 2 " *M 2 ±AM 2 ) (3 ' 18) 



t 2 =0 g^r|a=ao,X 2 =0 \ / \ ( SM\ ± AM; 

>,X 2 = 

which can be easily solved 



5a A / ^rU=a ,x 2 =o §|xU=a , x 2 = o \ / 5M\ ± AMi 
X 2 J " ^ ^U= QOlX 2 = o |f L=« ,x 2 =o J V ^ 2 ± AM 2 



(3.19) 



From this we can directly read off the allowed values of % 2 . 

For the cases of interest to us it is typically sufficient to determine the derivatives dMi/da, 
dMi/dx 2 to lowest order in a and x 2 • F° r example, in the previous subsection we have used first 
order perturbation theory, and determined the shift in energy to be 5E = x 2 /(my), accordingly 
dE/dx 2 = /(my). 

Let us confirm that this procedure corrects the behaviour of bounds for small m 2 ,. We 
consider a simple example where, in absence of hidden photons, both observables behave as 
simple power laws, ~ a 111 and ~ a" 2 , respectively, 

Mi = Mi\ x 2 =0 + 5 x2 Mi = Cl a n > + aa ni x 2 fiirn^) + 0( X 4 )- (3.20) 
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The term in 5 X 2 Mi is the hidden photon contribution calculated in first order perturbation 
theory (|3.ip using the potential (|2.3p . Note that for convenience we have written the correction 
with Cj and x 2 factored out, with the remaining factor represented by a function /j(my). 

From Eq. (j2.3|) we can see that 5V(r) term dies off at large my. Hence the function 
/i(my) —> as my — ?• oo. However, the perturbation actually grows towards smaller masses, 
and /i(my) tends to a constant limit. Before addition of the hidden photon we can write the 
Coulomb potential as 

V (r) = -— (3.21) 
r 

and after the hidden photon is added we can write 

V{r) = -^{l + X 2 ) (3.22) 

where we have taken the limit my — > 0. Hence the fine structure constant has essentially been 
redefined 

q^q(1 + X 2 )- (3.23) 

Therefore, for my —> 0, we have an alternative way to obtain the perturbation: we can 
simply insert the redefined a into the unperturbed expression. This yields, 

Mi = Ci a n * + n iCl a n >x 2 + 0(x 4 )- (3.24) 

Comparing (|3.20p with (|3.24p we can see that /j(my) — > Hi as my —> 0. 
Inserting (|3.20p into our general expression (|3.19p we find, 

i(5 A/2 n 2 SM 1 

^2 _ A/ 2 (a o ,0) A/i(a ,0) ^ ^ 

(ni/ 2 (my) - n 2 /i(my)) ' 

If both 5 Mi are consistent with we can obtain a bound on x 2 ■ To be conservative we 
simply add the moduli of the two individual errors, 

ni|AA/ 2 | , ra 2 |AA/i| 

X < IT — 77 v 77 777- 3.26 

|(ni/ 2 (my) - n 2 /i(my))| 

where |AM| is defined in (|3.4p . For transitions in atomic hydrogen and muonic hydrogen we 
use | AM* | defined in (|3"3j) . 

There are two interesting limit^l 

• my — > oo: Here /i(my), / 2 (my) — )• so the upper bound on x 2 increases and the bound 

(0-1-21) 

dies off. Indeed one finds that the functions /(m 7 ') decay as my at high masses 
where / is the lowest angular momentum value involved in the transition. A 2s — 2p 
transition would have I = and a 3p — 3d transition would have I = 1 etc. Our bounds 
effectively die off as 



,(i+0 

V 

Therefore, if we use transitions with higher values of / the bounds die off more quickly. 



X-m^r'. (3.27) 



my —7- 0: Here _fi(my ) — > m, / 2 ( m -7 / ) - > n 2 so the denominator tends to zero. Again the 
upper limit on x 2 increases and the bound dies off. 



2 We note, that in the regions where the bounds become very weak and \ = 0(0.1 — 1), perturbation theory 
breaks down. But these regions are typically excluded anyway. 
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Overall, we obtain the expected behaviour in the small and large mass limits. 

In Fig. [3] we show the correctly renormalised lsi/2 — 2s\j 2 (green) versus the naive bound 
(red). Here, and unless otherwise stated we will take M\ to be the 2s\i 2 — 8si/ 2 transition 
in atomic hydrogen, simply because it is experimentally measured to a high precision of 3 x 
1CP 11 eV [31]. This is similar to the theoretical error of 4.5 x KT 11 eV from and 5 x 10" 11 eV 
from the Lamb shift of the 2si/ 2 state. Added together with the AM(r p ) contribution we get 
an overall uncertainty of 2 x 10 eV at the 1 a level. 

The measurements themselves do not need to be from atomic spectra. The same technique 
works for any process which is affected by the hidden photon. 

Also note that this renormalisation procedure works trivially for same n transitions. If we 
consider just one measurement M, and note that the 0th order energy Mq vanishes, 

M = caV/ 2 (my) + 0(x 4 ). (3.28) 

M oc x 2 , so that 5a corrections obtained from a second measurement would produce negligible 
terms of 0(x 4 )- Therefore information from the second measurement is suppressed. This 
means that we can form properly renormalised bounds for same n transitions using only one 
measurement M and with only one error AM. 

4 New bounds 

4.1 Ordinary atomic hydrogen 

The bounds from 2si/ 2 — 2pi/ 2 (blue) and the properly renormalised lsi/ 2 — 2s\/2 bound (green) 
are shown in Fig. [3l We can see that neither of these penetrate new parameter space for hidden 
photons. However, they constitute the best pure tests of Coulomb's law in this region of 
parameter space. 

So far, all of our transitions have involved the lsi/ 2 or 2s 1 / 2 states. These have Lamb shifts 
with high theoretical uncertainty due to finite nuclear size effects. This is because the 0th order 
wavefunctions for s states are non zero at the origin and therefore penetrate the nucleus deeply, 
leading to large uncertainties from finite nuclear size effects. This will also be true for s states 
with n > 2. States with I > have 0th order wavefunctions which are zero at the origin and 
therefore have a small overlap with the nucleus. This means that the Lamb shift uncertainties 
are small. Therefore transitions between same n states with />0 (for example 3pi/ 2 — 3d>i/ 2 ) 
will have extremely small theoretical uncertainties. Unfortunately, these kinds of excited states 
are unstable and the experimental measurements have large uncertainties. We also note that 
the bounds in these transitions decay more quickly with energy (|3.27|) . causing them to die off 
well before they reach unknown parameter space. Therefore, barring significant technological 
advances, transitions between states with / > are not very useful for constraining x- 

4.2 Ions with Z > 1 

The Lamb shift bound for atomic hydrogen almost penetrates the unknown region, so we can 
look at hydrogenic ions with Z>\. The advantage here is that the characteristic energy of 
transitions scales as Z, so higher values of Z will move our bounds towards the right and towards 
the white region. The disadvantage is the higher theoretical uncertainties involved. This is due 
to an increase in the size of the nucleus, as well as a decrease in the Bohr radius (which scales 
as ^ and gives the characteristic length scale of the electron orbit). Unfortunately, this also 
causes the electron to penetrate the nucleus more deeply, which leads to greater theoretical 
uncertainties from finite nuclear size effects. Our bounds then move upwards and away from 
the white region. We can produce bounds for different values of Z to see which of these effects 
wins out. 
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Figure 4: \ bounds for the 2s\/2 — ?>Pi/2 transition in hydrogenic ions with Z = 2 (black), 
Z = 15 (green), and Z = 110 (yellow dashed, speculative). At the 2 a level the uncertainties 
are 3 x 10" 9 eV [21], 6 x 10" 4 eV J32JE3] and 8 eV [33] respectively. We can see that as Z 
increases our bounds move away from the white region and therefore become less useful. 



Ref. |24j gives experimental and theoretical errors for the helium ion 2si/2 — 2pi/2 transition 
of 0.16 MHz and 0.2 MHz respectively, giving AM = 1.5 x 10 -9 eV at the 1 a level. (Remember 
that for atoms without a proton nucleus, we do not have to consider large variations in nuclear 
size, and therefore use (|3.4p .) For Z = 15 the experimental [32] and theoretical errors [33] 
combine to give a value of AM = 3 x 10" 4 eV at the 1 sigma level. Finally we can go to the 
largest value of Z for which data is available. Ref. [33] gives a theoretical uncertainty for the 
Z = 110 transition of ~ 4 eV. As no experimental data exists we can only derive a speculative 
bound. Also, if we were to take into account hypothetical experimental data, then the errors 
would increase and the bound would weaken. The estimated sensitivity shown in yellow in 
Fig. [H is therefore an optimistic one, but it is still adequate to demonstrate the trend which we 
are trying to identify. 

We can see that as we go to higher values of Z, the increase in uncertainties cause our 
bounds to move up and away from the untested region, with our bounds actually becoming less 
and less stringent. We should expect this trend to be similar for other possible transitions, for 
example lsi/2 — 2si/ 2 . 

Note that due to the increase in the uncertainty of the Lamb shift in atomic hydrogen caused 
by the inconsistent values of the proton radius, the bound from hydrogen-like helium is actually 
slightly better. As soon as this inconsistency is clarified the trend identified above will probably 
hold again. 

4.3 Exotic atoms 

Exotic atoms have certain advantages over atomic hydrogen: 

• Pure QED systems like muonium and positronium may have smaller fractional theoretical 
uncertainties, as all experimental data is consistent with pointlike leptons [34] , This 
enables us to assume a pointlike effective nucleus, which eliminates the major source of 
theoretical uncertainty. 

• Most exotic atoms have larger reduced masses than atomic hydrogen, shifting our bounds 
to higher energies and towards the untested region. 

Disadvantages: 
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• Higher mass systems often have greater theoretical uncertainty (leptonic systems ex- 
cluded). This is because the larger reduced mass leads to a smaller Bohr radius of the 
system. Bound state orbits then penetrate the nucleus more deeply, which then leads to 
greater corrections due to finite nuclear size effects. 

• Similarly, hadronic atoms and atoms with Z>1 have larger nuclei than atomic hydrogen, 
again causing larger finite nuclear size effects. 

• Hadronic orbiting particles also interact with the nucleus via the strong interaction, 
which causes huge theoretical and experimental uncertainties. They do not produce strong 
bounds, but for completeness we briefly review them in App. [Cj 

4.4 Leptonic atoms 

4.4.1 Positronium 

We briefly note that positronium is not useful as the reduced mass is actually smaller than that 
of atomic hydrogen, and the uncertainties are much higher. For example the Is — 2s and 2s — 2p 
transitions are limited by large experimental uncertainties, which are caused by complications 
such as annihilation [35H37] . The uncertainties are around two orders larger, which gives a 
bound on x around one order of magnitude weaker. 

4.4.2 Muonium 

Ordinary muonium (fi + e~) and true muonium fj,~) are more interesting. 

Experimental results have already been produced for the lsi/ 2 — 2si/ 2 and 2s\/2 — 2>Pi/2 
transitions in ordinary muonium, but the resulting bounds do not improve on atomic hydrogen. 
We note that the reduced mass of ordinary muonium is almost the same (in fact slightly smaller) 
than atomic hydrogen, so we could only get better bounds if uncertainties were reduced. 

In fact the lsxn — 2s\m transition suffers from large theoretical errors associated with the 
effective Rydberg constant (|3.3p . This produces a fractional error larger than that of atomic 
hydrogen and consequently a weaker bound. 

The 2si/2 — 2pi/2 transition suffers from no such Oth order uncertainties, and indeed we 
expect the theoretical uncertainty to be smaller than atomic hydrogen since there are no finite 
nuclear size effects. However the experimental situation is not yet very good. Only lsx/2 states 
can be produced in large quantities, whereas the production of metastable 2su2 states is much 
lower [38]. This means that the 2s]y 2 — 2pi/2 transition has only been experimentally measured 
to a fractional precision of 1.5 x 10~ 2 at the 1 a level, leading to weak bounds. 

A recent article [39] suggests that true muonium can be produced and studied in the near 
future. This atomic system would be extremely useful, as it has a reduced mass of around 
2 orders of magnitude greater than atomic hydrogen, and we expect the theoretical errors to 
be low due to the absence of finite nuclear size effects. Since no experimental data has been 
produced yet, there has been no motivation for precise theoretical calculations. However, we 
can put together a rough theoretical estimate (see App. [B] for details) to form a speculative 
bound. This is shown as the black line in Fig. [5] This is encouraging as it penetrates new 
parameter space. However, one still needs to obtain a coherent experimental result, and hope 
that the experimental error is not so large that it causes the bound to significantly deteriorate. 

4.5 Muonic atoms 

The reduced mass of these systems is around 200 times larger than atomic hydrogen, shifting 
our bounds towards larger masses. 
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Figure 5: The black dashed curve shows a speculative bound on x f° r the 2s\i2 ~ %Pi/2 transition 
in true muonium. This is formed using only an estimate of the theoretical uncertainty of ~ 0.1 
GHz (see appendix [B]). If the experimental result is measured to a similar precision and agrees 
with theory, we will be able to form a bound similar to this and penetrate new parameter 
space. The green solid curve shows an actual bound formed from the 2sj^ 1 — 2p^ 2 transition 
in muonic hydrogen. There is a large theoretical error AM* = 1.2 x 10 -3 eV at the 2 a level. 
This is caused by the wide range of possible values of r p . The red dashed curve uses the muonic 
hydrogen transition to form a speculative bound. The error here is taken to be just the 2 a 
experimental uncertainty of 6 x 10~ 6 eV [25]. Finally, the blue curve shows a bound obtained 
by combining the measurement of the Lamb shift in ordinary hydrogen and the 2s^ 1 — 2p|^ 2 
transition in muonic hydrogen. Using these two measurements we do not need an additional 
determination of the proton radius. 



The 2sf / = 1 - 2pf,= 2 difference in muonic hydrogen is calculated to be 



E = — 



209.9779(49) - 5.2262 + 0.0346 



meV. (4.1) 



If we substitute in the most precise current value of r p = 0.8768(69) fm which is obtained 
from atomic spectra [23J, we obtain a theoretical value of E = —205.984(062) meV. The theo- 
retical uncertainty alone is quite high. Moreover, the theoretical value also deviates from the 
recently measured experimental value of —206.295000(3) meV |25| by around 5<r. This large 
discrepancy is bad for producing bounds, but it could be taken as a potential signal for new 
physics. We will consider this in the final Sect. [5j However, we can still form a bound if 
we use the (inflated) error given in (|3.5p . The r p variation increases the uncertainty to give 
AM* = 6 x 10" 4 eV at the 1 a level. The solid green curve shows the corresponding bound at 
the 2 a level. In the large mass region this bound is of comparable strength to the one from 
the Lamb shift in ordinary hydrogen. 

However if we just use the experimental uncertainty we can form a speculative bound Fig. [SJ 
This bound covers a similar region to the speculative bound obtained from true muonium (black 
curve) and penetrates the unknown region. 

If an independent and sufficiently precise value of r p could be determined - consistent 
with the muonic hydrogen extraction - the speculative bound could be turned into a real 
one. This provides us with motivation for seeking more precise, independent determinations 
of the proton radius. However, as we have discussed in Sect. O in order to avoid degeneracies 
this measurement should preferably originate from measurements at relatively high momentum 
transfer. At the moment the only obvious candidate process is electron scattering, although an 
increase in the precision by an order of magnitude may be challenging. 
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5 Muonic hydrogen anomaly* 

As already mentioned above a recent measurement [25] of the Lamb shift in muonic hydrogen, 
or more precisely the 2s^ 1 — 2p|^ 2 transition, deviates by more than 5cr from theoretical 
calculations combined with atomic spectra measurements of the proton radius. It is tempting 
to speculate that this deviation is due to a hidden photon [36]. In this section we will briefly 
investigate if such an interpretation is possible. 

The first observation is that the addition of the hidden photon increases the binding energy 
of the s-state compared to the p-state. This makes the Lamb-shift in muonic hydrogen more 
negative, in line with the observed effect. 

Encouraged by this we would like to do a fit using the Lamb shift in ordinary and muonic 
hydrogen to fit the proton radius and the kinetic mixing parameter x 2 of the hidden photon. 
We can use the same strategy as outlined in Sect. [3] just including the proton radius as an 
additional parameteiH. The hidden photon contribution to the Lamb shift is already given in 
Eq. (|3.8[) . For muonic hydrogen we just have to replace the electron mass with the muon mass. 
The lowest order dependence of the Lamb shift on the proton radius is given in Eq. (|3.6p . 

Unfortunately, the required values for x 2 are smaller than zero and since x is a real parameter 
this rules out a hidden photon explanation. Why is this the case? From Eqs. (|3.6|) and (|3.8|) 
we see that from the perspective of the (n = 2) Lamb shifts a shift in the proton radius by Ar 2 
is equivalent to a non-vanishing x 2 f° r 

4 2 

,, , , d t^TTL / 

Ar l = ~6X 2 ° 7 , 4 - (5-1) 
p (l + a m 7 /) 4 

where a Q = l/(am ) is the Bohr radius of the orbiting particle. 

This means that if we have non- vanishing x 2 > we measure effectively measure a smaller 
proton radius in the Lamb shift measurement. Obviously both the ordinary and the muonic 
hydrogen measurements are affected in the same direction. Now, one can easily check that the 
effect is actually always bigger for larger a a . In other words if there is a non- vanishing x 2 > 
then the observed proton radius in the Lamb shift of ordinary hydrogen should be even smaller 
than the one observed in muonic hydrogen. This is exactly the opposite of what is observed in 
Ref. [25]. 

Finally, we can use the same two measurements to form a bound, independent of electron 
scattering determinations of the proton radius. However, in light of the fact that the two 
measurements are not consistent with each other we need to inflate our uncertainty similar to 
Eq. (|3.5|) . Taking M\ as the Lamb shift in atomic hydrogen and using Refs. I22EE] we hav(S 
5 Mi = 3 x 1CT 11 eV and AM* = 10~ 10 eV at the 1 a leveB Taking M 2 as the muonic hydrogen 
case, and using the experimental and theoretical values from 14.51 we get 8M2 = —3.11 x 10~ 4 
eV and AMJ = 6 x 10 -4 eV. The corresponding bound is shown as the solid blue line in Fig. [5j 



6 Conclusions 

Atomic spectra can provide a powerful and clean probe of Coulomb's law at atomic length scales. 
In this paper we have considered a variety of transitions in atomic hydrogen, hydrogenlike ions 

*We are deeply indebted to B. Batell and M. Pospelov for noting a sign error in a previous version of this 
paper which led to different results. In the following discussion of the effects of hidden photons on measurements 
of the proton radius we will use an argument similar to theirs. 

3 As we are dealing with Lamb shifts (i.e transitions between states of the same n), changes in a are a 
subdominant effect (see Sect. I3.LT)| , 

4 The 5Mi are calculated using the CODATA [23] mean value for the proton radius of r p — 0.8768 fm. 

5 Note that since in this section we explicitly consider variations in r p in the bound, we could in principle ignore 
any contributions from variations in r p when calculating the theoretical uncertainty. However, as mentioned 
before, we will inflate our error to account for the inconsistency of the two measurements. 
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and exotic atoms, probing a wide range of different length scales. Currently, the best bounds are 
obtained from the measurements of the Lamb shift in ordinary hydrogen, a combination of the 
lsi/2 — 2si/2 and the 2s\/2 — 8su2 transitions in atomic hydrogen, as well as the Lamb shift in 
hydrogen-like helium ions. These bounds are shown in Fig. [1] and can be taken more generally 
as a constraint on deviations from Coulomb's law of the form Eq. (|2.3p . This significantly 
increases the range of pure Coulomb's law tests. 

Future measurements of true muonium could significantly increase the tested area 

towards shorter length scales and smaller deviations. Similarly, a bound from the 2s^ 1 — 2p|^ 2 
transition in muonic hydrogen could penetrate new parameter space, provided that a suitable 
independent value of the proton radius is determined. This gives additional motivation to seek 
further high precision determinations of the proton radius. 

The deviation from Coulomb's law, Eq. (|2.3p . and the bounds we obtain can also be in- 
terpreted as a probe of massive hidden photons kinetically mixing with the ordinary photon. 
With currently available data we find no improvement over existing bounds. We note, however, 
that these bounds are especially clean and model independent as they do not depend on the 
stability of the hidden photon or the presence/absence of certain decay channels. Moreover, 
above mentioned future measurements also have good potential to probe new parameter space 
for hidden photons, giving them significant discovery potential for new physics. 

Finally, we have briefly investigated a recent measurement of the 2s^ 1 — 2p|^ 2 transition 
in muonic hydrogen which found a value inconsistent with theoretical calculations and previous 
measurements of the proton radius. We have considered the hidden photon as a mechanism 
to explain this discrepancy and found that the measurements of the Lamb shift in ordinary 
hydrogen are in conflict with this hypothesis. 
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A Bounds on minicharged particles from atomic spectra 



Here we briefly examine bounds on minicharged particles (MCPs) produced by using atomic 
spectra. The existing bounds for MCPs are presented in Fig. [6j 

As noted in [7] the vacuum polarization caused by minicharged particles (shown in Fig. [7|) 
causes a modification to the Coulomb potential 



V(r) = VboulombM + 5V(r) M CP 
where the additional term is the Uehling Potential, 

2ae 2 



SV(r) M CP = — 
r 



3tt 



exp(— qr) 



1 



4m? 



2 m, 



1 + 



2mi 



(A.l) 



(A.2) 



Here m e is the mass of the minicharged particle, and e the charge in units of e. In the limits of 
high and low mass 5V(r) reduces as follows; 



SV(r) M CP 
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Figure 6: Existing bounds for MCPs (see, e.g. [H])- Note that the Lamb shift bound is produced 
using the method from [6]. We modify this method to consider other atomic spectra. 




Figure 7: Vacuum polarization diagram leading to a deviation from Coulomb's law in the 
presence of minicharged particles. 
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Figure 8: The lsi/ 2 — 2sij 2 (black) and 2s\i 2 — 2pi/ 2 [6] (green) bounds, both at the 2 cr level. 
The lsi/2 — 2si/2 transition is renormalised using 2s\/ 2 — 8s]y 2 as in the hidden photon case. 
For comparison we show previous bounds from purely laboratory experiments in orange. Our 
bounds do not penetrate new parameter space. 



and where 7 is Euler's constant. 

The potential decays exponentially in the high mass limit, as in the hidden photon case. 
However 5V{r)MCP does not die off in the low mass limit. This makes sense physically, as 
a massless MCP will still be distinguishable from the photon and therefore will constitute a 
non-trivial modification to the photon propagator, which will result in observable effects. In 
fact we find that the strongest effects are observed in the low mass limit. This can be seen from 
saturation of the bounds in Fig. [6J Hence we obtain strong bounds for 

m e < — (A.4) 
'0 

where Iq is the length scale of the experiment. 

From IA.3I we can see that the potential actually diverges logarithmically. We apply the 
same renormalisation procedure as we did with the hidden photon in Sect. 13.1.21 to get 



e < 



ni|AA/ 2 1 n 2 |AMi 
Ah + Mi 



\ nif 2 (m t ) - n 2 fi(m e )' 



Note the following limits; 

• m t — > 00: Here /i(m e ), f2{m t ) — > so e increases and the bound dies off. 

• m e — > 0: Here fi(m e ) and f2{m e ) both diverge logarithmically, but the divergences cancel 
and the denominator approaches a finite value as m e — > 0, giving us a saturated bound. 
This shows that, as expected, effects of the MCP are maximised in the low mass limit. 

We now apply this method to different atomic transitions. The analysis here is much briefer 
than in the hidden photon case, as we only consider transitions in atomic hydrogen. (The bounds 
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from systems with higher Z and exotic atoms are considerably weaker.) As noted before, higher 
excited states (i.e states which do not involve lsi/2 or 2s\j 2 ) are difficult to measure and have 
high experimental uncertainties. Therefore we only consider transitions involving the lsi/2 and 
2s]/2 states. In Fig. [8] we plot the Lamb shift 2s\i 2 — 2pi/ 2 EL and also a renormalised bound 
using the lsi / 2 ~~ % s i/2 and 2si/ 2 — 8si/ 2 transitions. The errors are the same as quoted in 
section 1370 We see that no improvement is found over existing bounds. 

B Estimate for the theoretical uncertainty in the 2si/2 — 2pi/2 
transition in true muonium 

The leading order contribution to the Lamb shift is proportional to the mass of the orbiting 
particle m . Therefore to get an approximate value of this transition we scale the atomic 
hydrogen value of by a factor ~ ^ to get ~ 1GHz. 

We note that the major part of the uncertainty in the atomic hydrogen case comes from finite 
nuclear size effects, and that these are absent in true muonium. We can therefore get a naive 
estimate of the uncertainty in true muonium by subtracting the finite nuclear size contribution 
from the atomic hydrogen value and scaling it up by ^ to get ~ 200 kHz. 

However this naive estimate is inadequate. The reduced mass of the system is now larger, 
so that hadronic and muonic vacuum polarization contributions are now much more important. 
These effects must receive more careful treatment. To leading order [47] 

Evpxf^Ymo (B.l) 
\m e J 

so that we need to scale up vacuum polarization contributions in atomic hydrogen by a factor 
(^■) 3 which gives us a much larger contribution of ~ 0.1 GHz. We therefore take this as our 
theoretical uncertainty, and set AM ~ 0.1 GHz. 

C Bounds from hadronic atoms 

We review this option briefly, and argue that hadronic atoms currently do not produce strong 
bounds. 

The existing candidates involve the tt~ , K~ , p + , and K~ particles orbiting a proton nucleus. 
In each case we find significant experimental uncertainties which immediately destroy bounds. 
For example with pionic hydrogen the experimental strong interaction shift of the lsi/2 ground 
state 08], 

e u = (7.120 ±0.008 ±0.007) eV (C.l) 

where the errors are systematic and statistical respectively. This essentially means that, for 
example, the 3p\/ 2 — lsi/2 transition would have an error of at least 10 -2 eV, which does not 
produce a useful bound. 

All transitions in kaonic and sigmaonic hydrogen have large uncertainties due to the deter- 
mination of the particle masses alone [49] . 

Transitions in antiprotonic hydrogen have large uncertainties due to both strong interaction 
shifts and annihilation |50j . 

This is before we take into account theoretical uncertainties, for example finite size effects 
from both nucleus and orbiting particle, which we also expect to be large. 

Overall we conclude that any benefits from the larger reduced masses of hadronic atoms are 
washed out by experimental uncertainties and QCD effects. 

6 Note that this has already been considered in Ref. [B] 

7 For values of the Lamb shift contributions in atomic hydrogen see, e.g. Ref. [30] . 
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